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As separable states are a convex combination of product states, the geome- 



Vidal and Tarrach are extended. Furthermore, it is proven that states in the 
set tangent to S at the maximally mixed state are separable; the set normal 
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constains, among others, all maximally entangled states. A canonical decom- 
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the closest product state to the maximally entangled state is the maximally 



mixed state. An algorithm is provided to find the closest product state. 
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PACS numbers: 03.67,03.65.Bz 

I. INTRODUCTION 

Entanglement has always been one of the quantum phenomena which sparked debates 
about the completeness and interpretations of quantum mechanics JT|, @, ||. It is funda- 
mental to teleportation [|], [[|, secure key distribution, |J, |7|], dense coding || and other 
applications. 

Obviously, therefore, properties of entangled states, including ways to determine if a state 
is entangled, are important. With regard to the latter, several authors have given sufficient, 
necessary, and necessary and sufficient conditions for a state to be entangled [|J, H, [JlO 
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IjJ . With the exception of Peres's partial transpose condition - if any partial transpose of 
a state has a negative eigenvalue, then the state is entangled - these conditions tend to be 
either impossible or hard to implement in all but low dimensional, low rank cases. This, of 
course, is an important starting point for the investigation of entanglement. 

In this paper we take a different approach. Rather than trying to determine if a state 
is entangled or not - the latter case being known as separable - we study the whole set of 
separable states and the whole set of entangled ones. Our approach is geometric and starts 
with the manifold of product states, which we denote by E. One reason for starting with 
E is that it is easy to decide if a density matrix is a product state by using partial trace. 
Another reason is that the totally mixed state ■hi is a product state. Lastly, E is of interest 
in itself. It is not convex or linear. But, as we shall show, it does have a nice geometric 
property: if a straight line intersects E in more than two states, then every point on the line 
is a product state. Sanpera, Vidal and Tarrach |12[proved this in the case of two qubits. 



We show that it is true in general. 

After establishing this geometric result, we turn to examining the set of entangled and 
the set of separable states. Our approach is to start at the totally mixed state, j^I and go 
in various directions. For instance, we show that all states which are in the set tangent to 
E at hi are separable. On the other hand the maximally entangled states are in the set 
which is normal to E at j^I. 

One of our main results is the fact that in the case of C n (g> C" the product state which 
is closest to a maximally entangled state is -kl, the totally mixed state. This surprising 
result, which we shall, unfortunately, show is false for systems with more than two particles, 
has an important consequence: there are no product states inside the ball of radius J 
centered at a maximally entangled state. 

It may seem to the reader that it would be hard to determine if a state is in the sets 
tangent or normal to E at -j^I and so the results just mentioned may seem to be hard to use. 
However, this is not the case. In fact it is easy to determine if a state is in the tangent or 
normal set. One uses a canonical, orthogonal decomposition of tq (N), the vector space of 



trace 0, Hermitian, N x N matrices, and the fact that every density matrix, Q, is uniquely 
expressible as Q = -^1 + H, where H G r (N). 

This fact and the canonical decomposition of r (N) are fairly standard mathematical 
fare; perhaps because itq (A) is the Lie Algebra su(N). However, it might not be well-known 
in the Quantum Information Community. Therefore, we shall present it in detail for general 
multi-particle systems. To get an idea of what is involved, consider the case of a density 
matrix Q on C"<g)C m . It has the orthogonal decomposition Q = j^I+H 1 ®^I+^I®H 2 +H 3 , 



where H\ G r (n), H 2 G r (m), and H3 G r (n) (g) r (m) . It follows from results in |L3| that 
Q is in the tangent space of S at j^I if and only if H3 = 0. Notice, it follows from 
Theorem 3 in the present paper that if Q is entangled, then H3 is not zero. We see from 
this decomposition, therefore, conditions under which entanglement occurs. As for Q being 
in the normal space of S at j^I, that happens if and only if Hi = Hi — 0. This follows, 
easily, from the orthogonality of the decomposition. To actually find the decomposition for 
a given state, one needs only to choose orthonormal bases for r (n) and r (m). The set of 
their tensor products then forms an orthonormal basis for r (n) ® r (m) . Given these three 
orthonormal bases, one then just takes inner products to get the decomposition. 

We finish this paper with an algorithm which appears, in the bipartite case, to give the 
product state which is closest to a given one. 

II. PRODUCT STATES 

We first establish some useful notation: 

T k (A) = {A I A is a Hermitian operator on C N and Tr(A) = /e} , (1) 

and 

VM{N) = {A I A is a density matrix on C^} , (2) 

where Tr denotes trace. Evidently, T>Ai(N) is the subset of T\ (A) consisting of positive, 
semi-definite operators. For iV = n\ • • • n p , we need to describe the set of product state 



density matrices on = C ni (g) • • • ® C np . Using the notation from JIB] , we have the 
product operator 

/i : nfri (m) -> n (I^n*) = n (iV) , given by (3) 
. . . ,Ap) = A X ®---®A P . 



The image of fi in n(iV), denoted by E, is a closed, embedded submanifold | 13[| . In 



particular, every point in £ has a tangent space and normal space. We shall use this later, 



but first we extend the result of Sanpera, Tarrach and Vidal [|12| about lines. Before doing 
so, we must point out that the set of product density matrices, denoted S + , in VM. (N) is 
just the image of [i restricted to iTj'P.M (n%) and so is a subset of S. Thus the next theorem 
holds a fortiori for product density matrices. 

Theorem 1 If A\ ® • • • ® A p — A and B G T\ (N) , then one of the following is true about 
the line r(t) = t (At ® ■ • • ® A p ) + (1 - t) B: 

1. It intersects S only at Ai ® • • • <g) A p , 

2. It intersects S at A\ ® ■ ■ ■ ® A p and exactly one other point, 

3. The line lies in S. In particular every density matrix on the line is a product state. 
Proof: 

Suppose B 7^ A is also a product, thus B = Bi<S> ■ ■ ■ <E> B p , and suppose t Ai <g> • • • <g> A p + 
(1 — to)Bi <S> ■ ■ • <8> Bp = Ci ® • • • ® C p for some t ^ or 1. In other words, suppose there are 
three products on the line. We need to show r{t) = tAi <g> • • • <8> A p + (1 — t)B\ ® ■ ■ ■ ® B p is a 
product for all t. If A 2 <g> • • • ® A p = B 2 <8> • • • ® 5 P , then we have r (t) = (tA x + (1 - ® 
^2 ® • ■ ■ <8> A p and so are done. Therefore assume A 2 <8> ■ ■ • <S> A p ^ B 2 <8> ■ ■ ■ <8> -Bp. Note, since 
Tr(A 2 ® • • ■ <g> A p ) = 1 = Ti(B 2 (g • • • ® we have that A 2 ® • • ■ <g> A p 7^ B 2 <g> ■ • ■ <g> 5 P is 
equivalent to them being linearly independent. 

Let {Ei} be a basis for t\ [rii). Then, using the Einstein summation convention of 
summing over repeated upper and lower indices, we have 



t A\Ei <g> A 2 ® ■ ■ ■ ® A p + (1 - to)-Bj^ ®B 2 ®---®B P = C[Ei ®C 2 ®---®C P . 

Since the i?j are linearly independent, this means that for each i 

t A\A 2 ®---®A p + (l- t )B\B 2 ®---®B p = C\C 2 ® ■ ■ ■ ® C p . 

If C % i = for some io, then = = B\°. For otherwise we would have toA l °A 2 ®---® 
A p + (1 — t )Bl°B 2 ®---®B p = 0, which would contradict our assumption that A 2 ® ■ ■ ■ ® A p 
and B 2 ® ■ ■ ■ ® Bp are linearly independent. Therefore assume CJ 1 ^ in which case we get 
t A ix (1 — t )B il 

-y-f A 2 ®---®A p + ± — °/ 1 B 2 ®---®B p = C 2 ®---®Cp. 

It follows from this and the linear independence of A 2 ® ■ ■ ■ ® A p and B 2 ® ■ ■ ■ ® B p that 
= t^L an d (i-y 1 = for eyery - i; - 2 guch that C n and c h are not Q 

Since we already know that A\ = B\ = if C\ = 0, we may conclude for all i that = 
and = 3Lci, where Cj 1 ^ 0. Thus A x = = Ad and ^ = = pd- 

Since Tr(Ai) = Tr(5i) = Tr(Ci) = 1, we have A = /i = 1 and so A\ = B\ = C\. This means 
we may write 



t A 1 ®---®A p + (l- t )B 1 ® ■ ■ ■ ® B p = d ® ■ ■ ■ ® (7, 



p 



as 



Ai ® (t A 2 ®---®A p + (l- t)B 2 ®---®B p )= A 1 ®C 2 ®---®C P . 

By induction we have that r(t) = tA\ ® ■ ■ ■ ® A p + (1 — t)B\ ® ■ ■ ■ ® B p is a product for all t. 
Hence if there are three matrices on the line r(t) which are products, then all are products. 

This actually proves the theorem, though one wants to know if the other two cases (one 
product or two products) can occur. In the next theorem we show that there are lines with 
only one product matrix on them. As for there being lines with only two, if that did not 
occur then every two elements in S would be connected by a line in S. This would mean E is 
convex. In particular every separable state would be a product state. However, this cannot 
be, for it was shown in that S + has measure in VA4(N); but the set of separable 



states has an open interior and so is not measure 0. Thus there are lines with exactly two 
product states on them q 

Before proceeding, we need to say a few words about the set of vectors tangent to E. 
Since E is the image of the embedding n, it follows that the tangent space at A — Ai<g>- • -®A P 
is o?/x(T'(A 1 ,...Ap)), where T^^.Ap) is the tangent space of n^ri(nj) at (A±, ...A p ) and dfu, is the 
Jacobian of /i. Thus the tangent space of E at A = Ai <g> • • • <g> A p is the vector space spanned 
by 

H 1 <g> A 2 <g> • • • <g> A p , A x <g> H 2 ® A 3 <g> • • • <g> A p , 
Ai <g> • • • <g> A p _i ® H p : Hi e t (m) 

Note, we distinguish between the tangent space of E at A and the set tangent to E at A. 
The former is the vector space of vectors which are tangent to E at j^I and so provide the 
directions which are tangent to E at A. The latter is the set of states which are obtained 
by starting at A and going in a tangential direction. Thus if T4 denotes the tangent space 
of E at A, then the set tangent to E at A is A + Ta- Similarly, the normal space of E at A 
is the vector space of vectors normal to E at A and the set normal to E at A is the set of 
states which are obtained by starting at A and going in a normal direction. If Ma denotes 
the normal space, then the set normal to E at A is A + Ma- 

Theorem 2 Let A = A\ <g> • • • <g> A p be such that each Ai = —I for i ^ i . If the line 
r(t) = tR + (1 — t)A is orthogonal to the tangent space o/E at A, then A is the only product 
matrix on the line. 

Proof: 

Suppose Q — Qi <g) • • • <S> Q p is also on the line. Since the line and tangent space are 
orthogonal, we have for all j and all Hj e r (nj) 

= (Q - A, Ai <g> • • • <g> Hj <g> • • • A p ) = (5) 
(Q 1 , A,)--- (Qj, H 3 ) ■ ■ ■ {Qp, Ap) - {A,, A,) ■ ■ ■ (A,-, H 3 ) ■ ■ ■ (A p , A p ) . 
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Now note that for % ^ i we have (Qi,Ai) = (Qi, ^-i) = ^-TrQj = and (A, A) = 



T (?7.j ). This can only happen if Q io = A io . However, if j ^ i then (|5|) reduces to 
(Q io ,A io ) (Qj,Hj) = for all Hj G r (nj), since in this case Aj = ^-1 and (^I,Hjj = 0. 
Thus (Qj, Hj) = for all Hj and so = — A? for all j ^ iq. □ 
We finish this section by showing every density matrix in the set tangent to X at the 
totally mixed state is separable. 

Theorem 3 Suppose Q is a density matrix and Q is in the set tangent to S at -j^I, then Q 
is separable. 



For Q to be in the set tangent to E at j^I, Q must be expressible as 

Q = — I + Hi <S> — I ® ■ ■ ■ ® — I H / ® H 2 ® • • • ® — I H (6) 

iV n 2 n p ni n p 

+— I® ■ .. ® _J_/(g)ff. 

If | ^j) is an eigenvector of ifj with eigenvalue Aj, then it is easy to see \tpi - • • ip p ) is an 
eigenvector of Q with eigenvalue -^(1 + n^Xi), where, as above, N — ni • • ■ n p . Since each 
Hi has trace equal to 0, the minimum eigenvalue of each Hi is 0, in which case Hi = 0, or is 
negative. Let /Zj be the minimum eigenvalue of Hi for each i. Then fi = -^(1 + J^i^^i) * s 
the minimum eigenvalue of Q. Since Q is a density matrix it is positive semi-definite and so 
< fi. This in turn implies j^I + —I (g) • • • <g) Hi (g) • • • ® — J is positive semi-definite and so a 
density matrix. In fact, if Ui ^ 0, then jjl + —I <S> ■ ■ ■ <8> Fj ® • • ■ <g> — J, where Fj = —t— x H^ 

iV Til Tip 

is a density matrix, since the most negative eigenvalue of Fi is ^ , = — — . 



Noting that £l + -^J ® • • • ® F ® • • • ® = ±1 ® ■ ■ ■ <S> (±1 + F;) ® • • • ® ±1, we see 

iV /fc^ / fcp ft^ v / ' / fcp 




Thus when j = i , (Eft becomes (Q io ,H io ) = (A io ,H io ) for all i/ G 



Proof: 



that Q is separable, for it is the convex combination 



P 1 P 1 1 

(1 - V rii l/Xil)— J + V n t \fx t \ (—1 ® • • • ® (— J + F») 




□ 
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III. A CANONICAL, ORTHOGONAL DECOMPOSITION FOR DENSITY 

MATRICES 



In this section we present a very useful orthogonal decomposition of trace 1, Hermitian 
matrices in terms of the totally mixed state and trace 0, Hermitian matrices. To begin, 
note that if Q is a trace 1, Hermitian matrix on then Q — -hi = S has trace and so 
is in the vector space r (N) . Also note that a trace matrix is orthogonal to -^J, since 
(jjI,S) = j^Tr(S) = 0. Thus we need to decompose tq(N), where N = n\---n p , into 
orthogonal subspaces. 

Let Z 2 = {0, 1} . The product of p-copies of Z 2 is denoted Z 2 . It is the set of all strings 
of length p of zeros and ones. For a = (a±, . . . a p ) G Z 2 take \a\ = Y7i a i- Thus |a| is the 
number of ones in a. For a G Zg, take u(a) to be the vector space u(a) = vi(a)®- ■ -®v p (a) , 
where Uj(a) = ^-1 if atj = and Vj(oc) = T (rij) if atj = 1. With these notational conventions 
we have, 

r (iV) = ©Li© H=MeZ ^(«). 

That this is an orthogonal decomposition of tq (N) follows from the fact that if A = A\ (g) 
• • • ® A p and B = B 1 ®- ■ ■ ® B pi then (A, B) = [T| = i (A h B { ) and the fact that H) = 
if H is trace 0. 

As can be seen from @) the subspace with k = 1, forms the directions which are tangential 
to E at jjl. Thus Q is in the set normal to S at j^I if Q = hi + S and S G © 2 ®\a\=h v(oi). 

To close this section we mention that the only role trace 1 plays in this decomposition 
is to determine the scalar multiple of I. In particular, any observable can be similarly 
decomposed. 



S 



IV. STATES THAT HAVE THE TOTALLY MIXED ONE AS THE BEST 
PRODUCT STATE APPROXIMATION 



Suppose Q G Ti (iV), where N = rii ■ ■ ■ n p . If P is the closest point on £ to Q, then the 
line joining Q to P is perpendicular to the tangent space of £ at P. Thus the only density 
matrices which could have the totally mixed state as best product state approximation (i.e. 
have jjl as the closest product state) lie in the set normal to £ at j^I. We showed in the 
last section how to characterize this set and shortly we shall show the maximally entangled 
states lie in it. But first we extend Theorem 3 by using the decomposition in the last section. 

Set T = ®\ a \ =1 v(a) and M = ® p k=2 ®\ a \ =k u(a) . 

Theorem 4 Every density matrix, Q, on C ni ® • • • <g) C np is uniquely expressible as Q = 
jfl + Qt + Qmi where Qt G T and G Af . In fact this is an orthogonal decomposition 
of Q. If Q is entangled, then ^ 0. If the totally mixed state is the closest product state 
to Q, then Qr = 0. If = 0, then Q is in the set tangent to £ at j^I . If Qt = 0, then 
Q is in the set normal to £ at j^I. 

Suppose N = n p and | ipi) is an orthonormal basis for C n . Modulo, U{n) <S> ■■■ <8> 
U{n) action (i.e. local operations), the maximally entangled state associated with this 
orthonormal basis is the projection onto ^5^" | with | ipi...ipi) being the p-fold 

tensor product of | ip-i)- Using standard tensor analysis notation, this projection is 



For what follows we need to compute (E, B) for an arbitrary trace 1, Hermitian matrix, 
B. Using the decomposition in the previous section, we know that B is the orthogonal sum 
of matrices B(a), a G 1F 2 , where B(a) = Bi(a) ® • • • <g) B p (a), with Bj (a) = ^1 if otj = 
and Bj(a) G r (n) if aj = 1. Thus we only need to compute (E,B(a)), since (E,B) is the 
sum of such quantities. 








1/n if ii 



otherwise. 




3 pi 



(7) 
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Since B (a) is a product (S(a))^"j Sp = {Bi(a))*l ■ ■ ■ (B p (a))^. Again using the Einstein 
summation convention and using the fact E 1 ^? = unless i\ = ■ ■ -% p and ji = ■ ■ ■ j p , we 
obtain 

(E,B(a)) =Ef;2(B 1 (a))il---(B p (a))r p , (8) 
= E-^ia))! ■ ■ ■ (B p (a))i 

When \a\ < p at least one of the ctj = and so at least one of the Bj (a) = -I. For this 

I if % ^ j 

one we have (Bi(a))j = < . Thus in this case (^j) becomes 

(E, B (a)) =Et:\{B 1 {a))\---{B p {a)% (9) 

= -(B 1 (a)y r --(B p (a))l. 

n 

Theorem 5 The maximally entangled state E is in the set normal to S at j^I. 
Proof: 

What it means for E to be in the set normal to S at j^I is that the vector parallel to the 
line which connects -kl to E is perpendicular to any vector tangent to £ at j^I. We have 
seen above that these tangent vectors are of the form 5^i a i=i E (a) ■ Hence we need to show 
that (E — jfl, X}[a|=i E = 0. This follows because 

(e - h, B («) ) = E «^ 5 (°)> - a7 v> b (°)> ( 10 ) 

\ |q|=1 / |a|=l 

and all the terms on the right hand side of (|10D are zero. Indeed the term (I, B (a)) = 0, 
since B (a) is trace 0. And the (E, B (a)) are because \a\ = 1. When \a\ = 1, there is only 
one factor in 5 (a) which is not -I, and so (E,B (a)) reduces to ^(Bj (a))], with Bj (a) 
being the one factor which is not -I. But (Bj (a))} is the trace of Bj (a), which is 0. □ 
Because of this theorem, we know there is the possibility that the closest product state 
to E is jrl. We shall now show that is in fact the case for two particles, but not the case 
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for more than two. First of all we note that the distance squared from E to any other state 
C is 



(E-C,E-C) = \\E\\ 2 - 2Re (E, C) + \\C\\ 2 • (11) 

1 1 1 1 1 2 1 

When C is the totally mixed state, j^I, this reduces to — j^. We need to compare 
this with the distance from £ to a product state. Suppose A is a product state. Then 
A = (~I + Rx) <S> ■ ■ ■ <8> (-/ + Rp) ,where Rj G r (n). Expanding this expression for A, we get 
A = jjl + Ylk=i ^2\ a \=k B(ct), where for each a G 1T 2 we have B(a) = Bi(a) ® • • • <S> -Bp(a), 
with Bj(a) = -I if ay = and Bj(a) = Rj if ay = 1. Recall from above that the B(a) are 

P i 

of them for each fc. Set -B^ = X]|a|=fc -^(°0- The 



mutually orthogonal and there are 

' k 

Bk are also mutually orthogonal. In terms of them, A = j^I + Y7k=i Bk- 
Using the mutual orthogonality of the Bk, we find 

(E,A) =/ E ,^l\ + J2(E,B k ), 

\ ' k=l 

k=l 

Similarly, ||A|| 2 = — + Ylk=i ll^fcll 2 - Hence it follows from ( [XT] ) that the distance squared 
from E to A is 

1 v p 

(E - A, E - A) = \\E\\ 2 - - - 2Re £ (E, B k ) + £ ||B fe || 2 . (12) 

k=l k=l 

1 1 2 1 

Recalling that the distance squared from E to the totally mixed state is ||_£7|| — jj, we 
see that the totally mixed state is the closest product state to E if and only if J^aUi ( II II 2 — 
2Re (E, Bk)) > for all choices of Ri, ...,R P G To (n) such that A is a density matrix. Let 
us first compute (E, Bk). For k < p we can use @ to get 

(£,fl fc ) = ^ (E,B(a)), (13) 

H=fc 



X)^i(«))l---(5p(«)): : - 



n 
H=fc 
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Suppose for a = (a±, ...,a p ) that acj = 1 for j — mi, ...,rrik. Then the summand in fll3|) 
becomes nP } k+1 {R mi )\ • • • {Rm k )\- Thus for k < p we have 

(E,B k ) = -^- k (R mi )l---(R mk )l. (14) 



n p+i- 

mi<---<m k 

For k = p we have (E, B p ) = l^HUY; ■ ■ ■ (i^ = J E^C^Oi " ' ' Wf- Hence 

2Re^(E,i? fc ) £ (ni? mi )^--(n J R m j: (15) 

fc=l fc=l mi<-<ra t 

Dispensing with the Einstein summation notation from here on out, we can rewrite (|Tj|) as 

(E,B k ) = (16) 
2 " * 2, 



EAR 1 + n(fli)j)) - 1)) + -Re 53(J2!)J • • • (R P )l 



i=l 1=1 i^j 

We also need to compute ||-£>a,-|| 2 . Fortunately, this is quite simple since the B(a) are 
orthogonal. In particular, we have for all k that 

\\B k f = \\B(a)\\ 2 =J2 \\Bi(a)\\ 2 ■ ■ ■ \\B p (a)\\ 2 , (17) 

|a|=fe |a|=fc 
— f ll^mJI 2 ---!!^ 



mi<— <mk 



Thus 



E™ 2 =^E E n||it! mi || 2 ...n||it: m J| 2 , (18) 

k=l k=l mi<—<mk 

= - p (f[(l + n\\R l \\ 2 )-l). 



i=i 



We can now rephrase our question as follows: is there a choice of Rj G T Q (n) such that 
A = (-1 + i?x) ■ • • (-J + i? p ) is a density matrix and 

-Allil + nWRtW 2 )-!) (19) 



2 2„ x • . ; 



(Edit 1 + - !)) - - Re E(^)' • • • to 



i=l i=l v£j 



12 



is negative. If so, then that product state is closer to the maximally entangled state than 
the totally mixed one. If not, then the totally mixed state is the closest product state to 
the maximally entangled one. 

Theorem 6 For a quantum system modelled on C n (g> C n the closest product state to a 
maximally entangled state is the totally mixed one. Because of this there are no product 



states within y ^r - , where N = n 2 , of a maximally entangled state. For quantum systems 
with more than two particles, the totally mixed state is not the closest product state to a 
maximally entangled state. 

Proof: 

Let us first consider the case of two particles. Thus p = 2 in (|I9|) and so it reduces to 

-^(l + n||i2i|| 2 )(l + npy 2 )--^ (20) 
2 n 
i=l 

n 

Expanding the first two terms, cancelling, regrouping and using the fact the R\ are 



Hermitian and so = {Rzjj we see that (P0|) equals 



\ (\\Ri\? + \\R2f ~ mi)\W\ + Re^W^jlj (21) 

+ ii J R 1 ii 2 ii J R 2 ii 2 -^(( J R 1 ): + (i? 2 )D- 

Since the Ri are trace 0, the last term in this expression is 0. On the other hand the first 
term is - — R2, R\ — R2), which is greater than or equal to zero, with equality only if 
Ri = i?2 = 0. Since the same is true for the middle term in (p0| ), we have proved that jrl is 
the product state closest to E in the bipartite case. 

Now suppose p > 3 and take all the Ri to be the matrix R with diagonal elements: 

l/2n Hi=j = l 
Ri = < -l/2n if i = j = 2 
otherwise. 
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When fll9f) has this R substituted into it, it becomes 



iff 1+ ±y + i_ 1 _2 (( ? )P+( i )P) 

up VV W n n 2 2 

For p = 3, this reduces to 8n 3 (— 12n 2 + 6n + 1), which is negative for n > 2. Since for 
such the quantity (1 + 4-) is less than |, for fixed n > 2, the expression (l + + - — 
1 — - ((|) p + (|) p ) decreases with increasing p and so is negative for all p > 3. 

In conclusion, the geometry of states was studied. Prior results by Sanpera, Vidal and 
Tarrach were extended and several new results were presented. These results are useful 
in the characterization of entangled states, which has been a difficult problem due to the 
complicated structure of the entangled space. 

Acknowledgement: The authors acknowledge the Office of Naval Research for the support 
of this work. 



V. APPENDIX 

In this Appendix an algorithm is presented to find the closest product state to an arbitary 
state. It was used, for example, to initially find that the closest product state to the 
maximally entangled state is the maximally mixed state. The algorithm converges very 
rapidly (often within a few iterations) for modest N, where N = n 2 . 

Given a density matrix C G C n <8> C n , consider the problem of finding 

arg min 1 1 C — A ® B 1 1 

A,B 11 11 

where 6 C" are density matrices. An outline of the algorithm is first provided, and 
later the details on the computation of A Q and B Q are given. 



A. Algorithm 

1. Set B — I n , . 
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2. Given B, find A Q as computed in Subsection B. 

• Set A = A + ±(l-Tr(A ))I n 

• Check if A is positive semi-definite (psd). If yes, then A satisfies a Lagrange 
Multiplier equation (see and discussion) skip to Step 3. If no, an A that 
satisfies a Lagrange Multiplier equation (see (|3~0D and discussion) can be found 
as follows: 

— Compute the eigenvalues of A D , represented by Ai > A2 • • • > A n and let 
$1, $2, • • • , $n be corresponding eigenvectors. 

— Set M = argmin m {m|(n — l)Aj + A m + l — Tr(A G ) > 0, j = 1, 2, • • • , n; j ^ m}. 

- For m = 1, 2, • ■ • , M, compute A m = A + ^(1 - Tr(A Q ) + \ m )I n - ^(1 - 
Tr(y4 ) + n\ m )$ m % n 

- Set A = argmin^ \ \C - A m <g> _B|| 2 

3. Given A find B Q . This is accomplished by considering the element transform T such 
that T(A®B) = B®A. Setting C = T(C), it is seen that \\C-A®B\\ = \\C-B®A\\. 
Then apply the methodology as in Steps 1. and 2. above with the appropriate changes 
to the transformed problem. 

4. Check \\C — A <g) B\\ for convergence. If no, go to Step 2. 

B. Solution for A Q and B 

Finding A D is equivalent to finding = 1, 2, ■ ■ • , n) such that 

aij = argmin ||Cy — ciijB\\. 

aij 

Now 

||Cjj — ciijB\\ 2 = Tr{(Cjj — aijB)(Cij — a^B)'} 

where ' denotes conjugate transpose. The above is minimized when 
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a 



a" 



" \\ B \? 

lTr^O + Tr^Q 



i?3, (22) 



" 2 ||fi|| 2 

Finding B Q can be easily accomplished with the same result as above, by considering the 
element transform T such that T(A <E> B) = B <g> A. Setting C = T(C), it is seen that 

\\C- A® B\\ = \\C-B®A\\, 

and the solution for b°- is by symmetry with (|22|) given by 

% ||A|| 2 

lIY(Q^)+Tr(AQ) 



6°. 
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2 PI, 

Remark The solution to A without the constraint Tr(A) = 1 for 

arq min 1 1 C — A ® B 1 1 
A 

where B, C are psd Hermitian (psdh) matrices, can be shown to be psdh. Hence both A Q and 
B are always psd. 

Consider now the solution to A with Tr(A) = 1 for 

arg min 1 1 C — A (8> B 1 1 

A 

where B, C are given density matrices. The Lagrange multiplier solution to 

Va {\\C - A® B\\ 2 + lTr(A)} = (23) 

can be shown to have the form A = A +^I n , where \j a is the gradient, Va=^- Substituting 
this form in the constraint equation, Tr(v4) = 1 gives 7 = 1 — Tr(A„). Since this is a convex 
minimization problem, an interior point solution to Lagrange multiplier gradient equations 
is guaranteed to be the minimum. If A is psd (A is Hermitian), then we will use this solution 
in the algorithm, although strictly speaking, it should be checked to see if it is an interior 
point. 
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The solution to A with Tr(yl) = 1 for 

arg min 1 1 C — A ® B \ \ 
Aepsdh 

where B, C are density matrices is not as straightforward if A = A a + (1 — Tr(A ))I n is not 
psd. In this case, we use the additional constraints x' Ax > and x'x = 1. The minimization 
problem is no longer convex. However, if the minimum is an interior point, then it necessarily 
satisfies the Lagrange multiplier gradient equations. As A = A Q + (1 — Tr(A ))I n is not psd 
but still minimizes | \C — A <g) B\ | 2 , then the solution lies on the boundary of the constraint 
x' Ax > 0, i.e. A is singular, hence x 1 Ax = 0. Now, the Lagrange multiplier problem becomes 
the solution to: 

Va {| \C - A ® B\ | 2 + 7iTr(A) + l2 x' Ax + ^x'x} = (24) 

V, -CI|C — A®5|| 2 + 7l Tr(A) +-f 2 x'Ax + j 3 x'x} = 0. (25) 

Eqn ( pop above merely shows that x is an eigenvector of A. Equation (^) imples that A 
has the following form: 

A = A Q + 7iJ n + 7 2 xx'. 
Now, 7i,72 must satisfy constraint equations: 

Tr(A) = 1 = Tr(v4 ) + + f 2 (26) 

and 

= Ax 

= A Q x + fix + 7 2 x 

= (A + (71 + i 2 )I n )x. (27) 

Equation (^) implies that x' is an eigenvector of A D . Denote the eigenvalues of A Q by 
Ai > A 2 • • • > A n and $i, $ 2 , • • • , $ n the corresponding eigenvectors. Let x = $ m , then 
Ax = implies 
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fi + 72 = -A m (A,,). (28) 

Using p7D and fl28|) we can show 

A = A + -^-r(l - Tr(y4 D ) + X m )I n - Ti(A ) + n\ m )$ m & m . (29) 

n — 1 n — 1 

We can show using ( |2"9"| ) that x'Ac = §' m A§ m = 0. It would appear, we have n possible 
candidates for the solution of A. each solution would be based on using a different <3> m in 
(p9|). However, we still must ensure that A is psd. The solutions for A are Hermitian. 

In order for A to be psd, we required that the eigenvalues of A (denoted by = 
1,2, • • ■ ,n) to satisfy jj > 0. Now, we can show if x = $ m then using (|29|) the following 
n — 1 conditions must be satisfied: 

(n - 1)7, = (n - l)Aj + A m + 1 - Tr(A ) > 0, for j = 1, 2, • • • , n, m. (30) 

Note Ai = maXj Xj. If the minimum occurs at an interior point then we are guaranteed that 
a solution for A must exist and 

(n - 1)A, + Ax + 1 - Tr(A D ) > (j - 1)A,- + A m + 1 - Tr(A Q ) 

for all j, m. Therefore 

(n - 1)A, + Ai + 1 - Tr(A ) > 

for all j. Thus we will successively check the largest eigenvalues of A Q and see which satisfy 
( p0|) (again Ai must satisfy (|30|) if an interior point minimum exists). The solution for A Q 
can be checked numerically to see if it is an interior point minimum. 
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